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Recent work has revealed that the wave function of a pure state can be measured directly and that 
complementary knowledge of a quantum system can be obtained simultaneously by weak measurements. 
However, the original scheme applies only to pure states, and it is not efficient because most of the data are 
discarded by post-selection. Here, we propose tomography schemes for pure states and for mixed states via 
weak measurements, and our schemes are more efficient because we do not discard any data. Furthermore, 
we demonstrate that any matrix element of a general state can be directly read from an appropriate weak 
measurement. The density matrix (with all of its elements) represents all that is directly accessible from a 
general measurement. 

For a group of blindfolded observers, who can each only touch one part of the elephant, to construct an 
accurate representation of an elephant, the group must combine information about different parts of the 
elephant (FIG la). In quantum mechanics, to have a complete description of a quantum system, in par- 
ticular, a complete knowledge of a quantum state, we must combine information about complementary aspects 
(FIG lb). 

An ideal (strong) quantum measurement of a certain observable only gives the probabilities of obtaining the 
eigenvalues of the observable, and the statistical results of the measurements reflect the diagonal terms of the 
density matrix in the eigenbasis of the observable. Measurements of different (complementary) observables 
provide the diagonal terms of the density state for different bases. If we perform measurements for a sufficient 
number of bases, we can reconstruct the density state, which is the idea of state tomography. 

An important difference exists between the tomography of a quantum state and its classical counterpart, e.g., 
the description of an elephant by blindfolded observers. In this classical example, the description (of the elephant) 
by each man is not exclusive, and the different aspects of perception can be simply combined to produce an overall 
description of the elephant (FIG la). However, quantum physics forbids simultaneous knowledge of comple- 
mentary observables because precise knowledge of a certain aspect necessarily implies uncertainty for the 
complementary aspects 1 . Therefore, one cannot simultaneously perform ideal measurements of complementary 
observables in quantum physics (FIG 2a). To determine the quantum state of a system, different measurement 
setups, each for a particular observable, are required. To determine the general state of a d- dimensional quantum 
system, at least d + 1 different experimental setups are needed. Each setup must be devoted to the measurement of 
one of the complementary observables. 

An alternative state tomography approach is possible. Instead of obtaining maximum information of a 
particular observable by an ideal measurement, we can perform a weak measurement 2 . The state of the quantum 
system is only slightly changed during its weak interaction with the measuring device; therefore, weak measure- 
ments of a set of complementary observables can be performed simultaneously 3 (FIG 2b). Thus, fewer experi- 
mental setups are needed for state tomography via weak measurements 4 . 

Weak measurements have been well described in the literature. Since they were first introduced by Aharonov, 
Albert and Vaidman 2 , weak measurements have been realised in experiments 5 " 7 , and have provided new insights 
into the study of paradoxes and fundamental problems in quantum theory 8 " 12 . Weak measurement has also been 
used as a practical tool for amplifying weak signals and signal-to-noise ratios 13 " 20 . Signal amplification via weak 
measurement usually involves a pre-selection and a post-selection that are nearly orthogonal, and the original 
formalism of weak measurement 2 ' 21 is not sufficient for the phenomena in this regime because it only retains the 
first-order terms of the interaction strength 22 " 24 . Extensions to the case of general pre-selection are given in 
(Refs. 25, 26). A framework that retains the high-order terms for general pre-selection and post-selection is given 
in (Ref. 27), and specific cases are further studied in (Refs. 28-30). A nonperturbative theory of weak pre- and 
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Figure 1 | (a) Compatible aspects are combined in classical physics, 
(b) while complementary aspects are combined in quantum physics. 

post-selected measurements is given in (Ref. 31). The general results 
for the weak measurement of a pair of complementary observables 
are provided in (Ref. 20). Related concepts, such as the contextual 
values 32 and modular values 33 , were also introduced. Weak measure- 
ments can also be used for state tomography. Lundeen et al. 4 directly 
measured the wave function of a pure state by a weak measurement 
followed by a post-selection. The scheme was extended to the mixed 
state by sequential weak measurements of pairs or triple products of 
complementary observables 34 . The results of the weak measurements 
could also be interpreted in terms of the Kirkwood-Dirac representa- 
tion via complex joint probabilities 35,36 . Reviews and additional refer- 
ences on weak measurements can be found in (Refs. 31, 37, 38). 

The results in (Ref. 4) demonstrate that the wave function of a pure 
state can be directly obtained from a single experimental setup, and 
that one can directly measure both the absolute values and the phases 
of the coefficients of a pure state in a certain basis. However, the 
scheme proposed in 4 applies only to pure states, and it is not efficient 
because most of the data are discarded due to the post-selection. Here, 
we propose tomography schemes for both pure states and mixed states 
via weak measurements. Our schemes are more efficient because we 
do not discard any data. We also show that any (diagonal or off- 
diagonal) element of the density state can be directly determined as 
the average pointer shift in an appropriate weak measurement. 

Results 

The idea. In the scheme proposed in (Ref. 4), the wave function to be 
measured is a continuous function of the position. The essence of the 
scheme is reviewed in finite dimensional Hilbert space as follows 34 ' 39 . 
Suppose a quantum system with a d- dimensional Hilbert space is in 
an unknown state 



d-i 



(i) 



where the coefficients in a certain basis {|a z )} need to be 
determined. The scheme for measuring the coefficients 
consists of a series of weak measurements of the observables A f = 
\cii) (fl/|, each followed by the same postselection, i.e., a projection 

onto the final state \b 0 ) = -j^yj'^aj). The weak value of the 
yd 

observable A,- is defined as 

(fc 0 |af)(fl,#) 1 



<&<#> Vd(boW 



(2) 



Both the real and the imaginary parts of the weak value W f have 
physical meanings and can be determined experimentally because 
they correspond to the average shifts of the pointer position and the 
momentum, respectively. From Eq. (2), we know that the weak values 
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Figure 2 | (a) A strong measurement only reveals a certain aspect of the 
quantum state, (b) while a weak measurement enables us to perceive 
(measure) incompatible aspects of the quantum state simultaneously. 

Wj are directly proportional to the coefficients ij/j that we want to 

measure and that the factor —= can be determined by the 

Vd{boW 

normalisation condition (up to an unphysical overall phase of \\//)). 
Therefore, we obtain a direct measurement of the coefficients thus 
a direct measurement of the wave function 

The essential point in this scheme is the choice of the post-selected 
state | b 0 ). Its overlap with each basis state | a t ) has the same magnitude 

and phase; therefore, the factor K = \jj i l ' W i = y f d(b 0 \\jj) does not 
depend on i and can be determined by normalisation. Because the 
change of the system state due to a weak interaction with the pointer 
is negligible, the success probability of post-selection is given by P = 
\(b 0 \\l/)\ 2 . Therefore, only a fraction P of the data is retained, and the 
remainder is discarded due to the post- selection. When the dimen- 
sion d is large, as in the case of a continuous wave function, the 
majority of the data is discarded. 

State tomography of a pure state. The part of the data that 
corresponds to the failure of the post-selection can be retained. 
We replace the final post- selection by a complete projective 
measurement onto the basis states {{bj}}. The inner products ftp = 
(bj\ai) are fixed when the two sets of the basis states are chosen. We 
organise the data according to the final state: if the final state is | bj), 
then the weak value of A f = \a t ) (a z | is given by 



(bj\ai){ai\^) _ (bj\ai) 



h 



•A, 



(3) 



Therefore, for each fixed j, the weak values of different A, give the 
relative ratios of the coefficients i//,: 



(bj\a,) h 



(4) 



Because {bj\\jj) does not depend on /, the coefficients i/^ are directly 

Un- 
proportional to — — , where the weak values Wj* are read- outs directly 

Pji 

from the experiments and the constants /? ; - f = (bj\a t ) are fixed when 

the two sets of basis states are chosen. The super index (j) in ij/^ 
denotes the coefficients ^ that we obtain from the jth subset of the 

data. Ideally, should not depend on j; however, they may depend 
on j due to imperfections and noise in the experiments. The 
magnitude of the factor (bj\i//) = e l(p ) can ^ e removed by 
normalisation: 
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The phase factor (pj has no physical meaning because it corresponds 
to the overall phase of the pure state we are attempting to measure. 
This factor can be removed by a convention. For example, we require 
xj/o > 0 (or the coefficient with the largest magnitude be positive). 
Using this approach, we retain all of the data and separate the data 
into d subsets according to the final states {\bj)}. For the jth subset of 
the data that corresponds to the fixed final states we have an 

estimation of the coefficients from Eq. (4) or (5). A different 
subset of the data gives a different estimation for the set of 
coefficients {i/^J. The differences between the estimations indicate 
the amount of error or noise in the experiment. 

When the basis {\bM is mutually unbiased to the basis {\a t )}, i.e., 

ftp = (bj\ai) = (5) can be simplified to 



(6) 



The original scheme proposed in 4 corresponds to the case in which 
| b 0 ) = ^ f=o I a i) ' an d only the data corresponding to the success- 
ful post-selection of \b 0 ) are retained. Here, we retain all of the data, 
and each subset of the weak values corresponding to a fixed final state 
gives an estimation of the state to be measured. The probability of 
obtaining the final state | bj) is given by |(fy|i/f)| 2 because the change of 
the system state during the weak interaction is negligible. When the 
state | to be measured is almost orthogonal to a certain final state 
\bj), then the relative frequency to obtain the final state \bj) in the 
post-selection is small. However, the corresponding weak value has a 
large magnitude, which could result in a large signal-to-noise ratio 
when the dominant noise is due to systematic error or imperfections. 

In the above scheme, weak measurements for a complete set of 
projective operators {A* = | a t )(cii\ } are needed. For the tomography of 
a pure-state wave function \\//), weak measurements of a single 
observable instead of a set of observables are sufficient. This single 
observable can be chosen as follows (an alternative choice is pre- 
sented in the Methods). We consider a pure state \(p) such that its 
overlap with each postselected state is nonzero, i.e., (bj\<p) ^ 0 for 
each;'. We perform a weak measurement of = \<p)(<p\. When the 
post- selected state is | bj), the weak value of P^ is defined as 

(bj\p,w) _ (bj\<p)(<p\M 



w,- 



Therefore, 



W, 
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Here, rjj = ^ is completely determined from the experimental 

data and our choice because Wj is directly obtained from the experi- 
ment and both \(p) and {| bj)} are fixed by our choice. From Eq. (8), we 
can reconstruct the state 

(9) 



W)=Y,Kn,\*})> 



where K is determined by the normalisation condition up to a phase. 

State tomography of a mixed state. To this point, we have only 
considered tomography schemes for a pure state. If the state we 
like to measure is a general mixed state p, we also have a 
tomography scheme via weak measurements. For this case, we 
perform weak measurements for a complete set of projective 
operators A z = |fl,-X<i/|, followed by a final projective measurement 
onto the basis {\bj}}. For the initial state p and the final post-selected 
state \bj), the average shifts of the pointer position and the 
momentum in a weak measurement of A { correspond to the real 
and imaginary parts of the weak value 25,27 



H(\bj){bj\)(\a 
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(10) 



From Eq. (10), one can express the matrix elements of p either by the 
basis {\a t )} 



or by the basis { | bj)} 



ik 



(11) 



(12) 



where the probability Pj = (bj\p\bj) corresponds to the frequency of 
obtaining the final state | bj) in the experiment. The weak values W),- 
and the probabilities Pj are directly accessible from the experiment. 
The relationships in Eqs. (11) and (12) show that the matrix elements 
of the mixed state can be written as linear summations of the weak 
values that are directly accessible from the experiment. 

Partial state tomography. Using the idea of weak measurements, we 
can also selectively measure some of the matrix elements (diagonal 
and off- diagonal terms) of a state directly. This is especially useful 
when we are only interested in one or a few matrix elements of the 
state and we do not need to perform state tomography for the entire 
density matrix. Suppose we are interested in measuring a particular 
matrix element (a\p\ b). We consider two different cases according to 
whether \a) and | b) are orthogonal. For the first case, suppose (b\a) 
0. Here, we can simply perform a weak measurement of \a){a\, with a 
post-selection onto the state \b). The weak value is defined as 



W~- 



(b\a)(a\p\b) 
(b\p\b) ' 



Thus, 



{a\p\b) = 



(bW) 
(b\a) 



W. 



(13) 



(14) 



Wis the weak value, and (b\p\ b) is the probability of success for post- 
selection. Both values are directly obtained from the experiment. 
{b | a) is only a fixed factor. Therefore, (a\p\b) is determined by 
Eq. (14). For the second case, suppose that (b\d) = 0. We then 
perform a weak measurement of the observable C = \c) (c\ with 
\c)= ^ (| a) + \ b)) and a follow-up projective measurement onto a 

basis including both \a) and \b) as the basis states. We must only 
retain the data when the final state is either \d) or | b). If the final state 
is | a), then the weak value of C is given by 



W-- 



(a\p\a) 



1 + 



(b\p\a) 
(a\p\a) 



If the final state is | b), then the weak value of C is given by 

„_ (b\c){c\p\b) _l( {a\p\b)\ 
(b\p\b) 2\ + (b\p\b)j- 



We have 



(b\p\a) = (a\p\a) (2W-1), 
(a\p\b) = (b\p\b) (2W-1). 



(15) 

(16) 

(17) 
(18) 



The post- selection probabilities (a\p\a) and (b\p\b) and the weak 
values W and W are directly accessible from the experiment; there- 
fore, we can determine the matrix element (a\p\b) from either 
Eq. (17) or (18), which should agree with each other. Any discrep- 
ancies between them could be used as indicators of the noise and 
imperfections in the experiment. 
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Experimental scheme. Here, we discuss how our schemes can be 
realised experimentally. For the weak measurement of an observable 
A, the interaction between the measuring device and the system is 
generally modelled by a Hamiltonian H=gA®pd(t — t 0 ), where p 
denotes the momentum of the pointer (measuring device). Because 
the interaction in a weak measurement does not change the state of 
the system significantly, weak measurements of several observables 
(whether they are commutative or not) can be performed 
simultaneously. The interaction of the simultaneous weak 
measurements of a set of observables {AJ can be introduced by the 
Hamiltonian: 

H=Y^gi A i®PAt-to). (19) 

i 

Here, pi denotes the momentum of the /th pointer that is coupled to the 
observable A,-. {A,-} is an arbitrary set of observables of the system. They 
could be a set of commuting observables (for example, A z = \a t ) (a,-|) or 
even a set of complementary observables. Suppose the initial state of 
the system is p s , and the initial state of the pointers is p d = p x ® P2 ® 
where p t denotes the initial state of the /th pointer. Suppose the 
post-selection of the system is a projection IT ; (for example, IT ; = | bp 
(bj\). Then the average position shift dq t and the average momentum 
shift dpi of the /th pointer are given by (see the derivation in the 
Methods) 

5q i =g i ReW ji (20) 



3p i = 2g i ImW ji (Ap i ) 



(21) 



The weak value Wp of A b for a general initial state p s and a general 
postselection ITj, is given by 20 

Tr(UjA iPs ) 



W jr - 



Tr{U jPs ) 



(22) 



The shifts of different pointers can be read out simultaneously, while 
either the position shift or the momentum shift is recorded at each 
time for each pointer. With the same Hamiltonian, we can obtain all 
of the weak values Wp. From these weak values, we can obtain all of 
the elements of the density matrix p s according to the state tomo- 
graphy strategy proposed in this article. 

Example of application. The state tomography schemes proposed 
here can also be used to detect a tiny parameter encoded in a 
quantum state. For example, to detect a tiny phase difference 3 
picked up by two orthogonal states of a qubit (a photon) that 
passes through a certain medium (a birefringent crystal), we could 

prepare a qubit in the initial state — = (|0) + |1)), and seek to 

V2 

determine the tiny phase 3 encoded in the outgoing state 

—j=i (|0) + e t3 \l)). To determine 3, we can perform a weak 
v 2 

measurement of the observable |1) (1| by introducing an 

interaction H =g\ 1) (1 1 ®p3(t — f 0 ) between the qubit system and a 
measuring device and a follow-up post-selection of the qubit system 

onto the state — ( 1 0) — 1 1 ) ) . The weak value is W = 1 — i — , and the 
v2 3 1 

average shift of the pointer momentum is given by dp = — 2g(Ap) — . 

Although the phase 3 is small, the average shift of the pointer 
momentum is not necessarily small, and this has certain 
advantages over the quantum-entanglement based scheme in 
Ramsey interferometry 40 , which relies on the use of GHZ states of 
many qubits, which are an expensive resource. 

Discussion 

In this article, we have proposed several schemes for state tomo- 
graphy via weak measurements. An advantage of our schemes is that 
we need to measure fewer observables compared with the standard 



method of state tomography based on ideal measurements. For 
example, a simple standard state tomography requires ideal measure- 
ments in at least d + 1 different bases; therefore, one needs to use at 
least d + 1 different experimental setups. Here, our tomography 
strategy requires weak interaction in one basis ({ | a t )}) and post-selec- 
tion in another basis ({ | bp}), and the same interaction Hamiltonian is 
used for all data collection. 

Compared with the scheme in (Ref. 4), our strategy applies not 
only to pure states but also to mixed states. Our strategy is more 
efficient because we retain all of the data and thus reduces the num- 
ber of experimental runs needed. Compared with the schemes in 
(Ref. 34) in which a density matrix is directly measured via sequential 
weak measurements of pairs or triple products of observables, our 
schemes are based on a single-time weak measurement followed 
by a strong complete projective measurement. Our interaction 
Hamiltonian in Eq. (19) is straightforward and differs significantly 
from those in the schemes of (Ref. 34). In the previous schemes 4,34 , 
the observables A z = {a t \ require separate weak measurements at 
each time for a particular i. In our schemes, the set of observables {A,-} 
can be measured simultaneously via the interaction introduced by 
the simple Hamiltonian in Eq. (19). 

Interestingly, any matrix element of a density state is directly 
accessible from a suitable weak measurement. Ideal measurements 
give probabilities and expectation values, which are only linear com- 
binations of the matrix elements of the density operator. Based on all 
of these facts, it is natural to draw the conclusion: what is and can be 
really measured in a general (ideal or weak) measurement is the 
density matrix (with all its elements), and only the density matrix! 

Methods 

An alternative scheme for the tomography of a pure state via the weak 
measurement of a single observable. In the main text, we see that it is sufficient to 
perform a weak measurement of a single observable for the tomography of a 
pure-state wave function This single observable could be chosen alternatively as 
follows. We introduce a single observable 



A=Y,k\a i )(a i \, 



(23) 



where the eigenvalues kj are not degenerated, i.e., k { ^ kj if i ¥= j. We simply perform a 
weak measurement of the observable A and retain all of the data corresponding to 
different final states \bj). When the final state is \bj), the weak value of A is given by 



Wj- 



0#> 



From Eq. (24) we have 



(24) 



(25) 



We define a unitary matrix /? with the matrix elements given by the coefficients ftp and 
introduce the diagonal matrices k = diag{k 0 , k l9 ' ' ' , k d -i\, w = 
diag{w 0 , Wi, ••• ,w d -i\ and the vector ijj = (\J/ Q ,\l/ l ,---,il/ d _ l ) T , where T denotes the 
transpose. Then, Eq. (25) can be written as 

(Pk-wP)\j/ = 0. (26) 

With the notation M = (Ik — wfi, we have Mxjj = 0. Thus, the state tomography in this 
case is to find the vector \Jj corresponding to the zero eigenvalue of M. The matrix M 
depends on w, which is determined by the weak values that are directly read out from 
the experiments. Due to noise and imperfections in the experiments, matrix M may 
not have a zero eigenvalue. Instead, multiplying M\]j = 0 by M 1 " from the left, we have 
(M f M)^ = 0. (27) 
Therefore, the vector \j/ is the eigenvector in the kernel space of the non-negative 
matrix M f M. The vector \Jj is uniquely determined if M*M has a one-dimensional 
kernel space. It is not determined, and another set of measurements must be 
performed if M f M has a kernel space of more than one dimension. However, it is also 
possible that the matrix M^M determined from the experiment may not have a kernel 
space. In such a case, we choose the most likely one, i.e., we choose the vector \j/ as the 
eigenvector corresponding to the smallest eigenvalue of M + M. The value of the 
smallest eigenvalue of M f M indicates the amount of noise and imperfection in the 
experiment. 
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Derivation of the position shift and the momentum shift. Here, we derive the 
position shift and the momentum shift in (20) and (21), respectively. The time 
evolution operator corresponding to the interaction Hamiltonian in (19) is given by 
U = e~ l ^2,i giAi ® pi . If the system is successfully post-selected by Tip which occurs with 
a probability (to the first order of g r ) of 

Pj = Tr((lL } ®i)Up s ®p d rt) 

( ^ A ( 28 ) 

= Tr(U jPs ) ll+lY^gilmWjiipi)] 

the final (unnormalised) state of the pointers is given by (to the first order of g { ) 

f/ d = Tr,((TI j ®l)Up t ®p d rt) (29) 




Here, Tr s (Tr) denotes the trace over the system (whole) Hilbert space. The average 
position shift of the z'th pointer, conditional on the successful post-selction by is 

given by Oqi = , - y — {C^ij.K straightforward calculation to the first order of 
ogives r v ) d) 

dq^giReWji+gilmWjidpi-iPi), (31) 

When the initial states of the pointers are well-behaved states, such as Gaussian states, 
the second term on the right-hand side of (31) vanishes, and we immediately have 
(20). Similarly the average momentum shift of the z'th pointer, conditional on the 

Tr (pip'd) 

successful post- selection by Ilj, is given by bp { = — , f , —{pi), which yields (2 1 ) to 

Tr \Pd) 

the first order of g t . 
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